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Elastic  Stresses  and  Deformations  Produced 
in  a  Semi-infinite  Elastic  Solid 
by  a  Point  Source  of  Heat  beneath  Its  Free  Surface 


For  the  purpose  of  solving  our  problem,  we  have,  first  of  all, 
deliberately  set  up  the  displacement  function  , 

0  —  LiL?  :*&.  (s/\g-  a|*  +  7-  >/  O  +  //)’  +  r‘)  (19) 

'  l  —  v  4xXo 

aiming  to  fulfill  the  Poisson  equation 

;  '  :  '"«(3^S-  +  £)V^-  uo), ;■ 


where  the  temperature  field  T,  aroused  by  the  point  source  of  heat, 

amounts  to .  : 

r  =  _2—  (— _ 1  — L_ -) 

.  iirli 'i/  \s  —  h\*  +  r‘  .-v  (v  i-  h):  +.  (17) 

; In  order  to .warrant  the  unloaded, conditions  of  the  free  surface  of 
■  the  semi-infinite  elastic  solid,  which,  in  general,  will  not  always  be. 
fulfilled  by ■ the . sought  displacement  function  0,  we  have  resorted  to 
a  search  for  a  second  displacement' function,  whose  derivation- was  carried- 
.  out  in  another  paper  by  the  author  entitled  "Elastic  stresses  and  deforma- 
,  tions  produced  in  a  semi-infinite  solid  by  the  application  of  tangential' 
traction  to' its  free  surface"  (LI-H5&EH  HS&EH-PAO/ACTA  MECHANIC*  SINICA, 
Vol..  3,  No.  2,  1959,  120).  “  • 

When  we  superpose  the  stresses  as  well  as  the  displacements  corres¬ 
ponding  to  the  two  above-mentioned  displacement-  functions,  we  obtain  the 


following  expressions  as  the  solution  to  this  thermal  stress,  problem; 
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All  the  above  expressions  for  the  stresses  and  the  displacements 
have  been  numerically  evaluated  with  appropriate  data  and  discussed  in. 
detail  with  accompanying  graphs  in  Figures  1  to  18  (original  text). 

The  variations  of  these  expressions  representing  the  complete  solution 
to  this  .thermal ':elastic  boundary  value  problem  are  shown  by  the  curves 
,in  Figures  19.  to  36- (original  text). 
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Similarity  Structure  of  Vorticity  Fluctuation 
and  the  Theory  of  Turbulence 

The  fluctuation  motion  in  a  turbulent  flow  is  considered  to  consist 
in  the  random  motion  of  vortices  governed  by  the  equations  of  turbulent 
vorticity  fluctuation  derived  by  taking  the  curl  of  the  equations  of 
velocity  fluctuation.  -In  the  downward-moving  system  at  mean  velocity 
of  the  fluid,  the  vortices  in  the  fluid  are  assumed  to  have  similarity 
structure.  In  other  words,  according  to  the  present  theory,  dynamically 
similar  solutions  are  chosen  for  the  equations  of  vorticity  fluctuation 
in  a  moving  system  of  reference. 

Special  problems  are  solved  for  turbulent  flows  limited  by  channel¬ 
like  boundaries,  along  a  semi- inf inite  plate  and  through  a  circular  pipe, 
and,  for  free  turbulence,  in  the  two-dimensional  and  axially  symmetrical 
wake3.  It  is  shown  that  by  neglecting  the  terms  for  the  degree  of  vis- 
.  cosity  in  the  equations  of  vorticity  fluctuation  in  the  moving  system, 
the  dynamical  similarity  condition  leads  to  von  Karman’s  well-known 
similarity  regulations  for  boundary-limited  flows  in  the  central  region 
of  the  wakes;  while  in  the  outer  region  within  the  wake  interior,  where 
the  effect  of  viscosity  cannot  be  ignored,  the  similarity  condition  yields 
the  experimentally  established  relation  that  the  Reynolds  apparent  shear¬ 
ing  stress  component  is  proportional  to  the  transversal  gradient,  of  the  . 
mean  flow  velocity; 


The  requirement  that  the  turbulent  vorticity  structure  should  have 
dynamical  similarity  leads  to  the  conclusion  that  all  the  components  of 
the  Reynolds  stress  are  mutually  proportional  to  each  other  in  the  flow 
region  where  a  definite  type  of  similarity  prevails.  This  theoretical 
result  has  experimental  verification  both  in  the  boundary-limited  flows 
a»ndi  in' the  two-dimensional  wake.  The  velocity  correlations  between  two 
distinct  points  should  also  be  independent  of  their  position  within  the 
flow  region  having  the  same  type  of  similarity,  so  long  as  the  relative 
displacement  between  the  two  points  remains  the  same.  This  also  agrees 
qualitatively  with  the  observations  available. 

The  application  of  the  condition  of  dynamical  similarity  gives  us 
a  length  A,  characterizing  the  size  of  vortices  which  compose  the  tur¬ 
bulent  fluctuation.  In  the  Karman  type  of  similarity,  A  depends  only 
upon  the  space  derivatives  of  the  meaii  velocity  and  is  independent  of 
the  fluid  viscosity,  while  A  of  the  other  type  of  vortices  in  the  outer 
part  of  the  wake  interior  is  dependent  upon  the  coefficient  of  viscosity. 
The  size  of  the  latter  kind  is  much  smaller  than  the  former.  These  two 
kinds  of  vortices  in  turbulent  motion  have  long  been  anticipated  in  theo¬ 
retical  investigations  of  turbulence. 

In  view  of  the  fact  that  the  condition  of  dynamical  similarity  can 
only  be  applied  within  a  definite  flow  region  and  that  not  only  in 
different  flow  problems  but  also  in  different  flow  regions  of  the  same 
problem  there  are  different  similarity  solutions  as  in  the  central  and 
outer  regions  of  the  wake  interior,  we  can  see  the  limited,  nature  of  the 
application  of  the  similarity  concept  in  the  turbulence  theory. 


Hence  more  rigorous  solutions  of  the  original  vorticity  fluctua¬ 
tion  equation  should  be  sought  in  order  to  connect  the  different  types 
of  similarity  solutions  present,  for  instance,  in  the  wake  problem. 

The  non-dimensional  equations  of  vorticity  fluctuation  are  still 
non-linear  after  the  similarity  transformation.  It  is  difficult  to  ob¬ 
tain  their  solutions  and  (based  upon  them)  to  build  the  velocity  corre¬ 
lation  functions  between  two  distinct  points  by  the  space  average 
method.  Since  Taylor's  scale  of  microturbulence  A  is  defined  by  a 
double  velocity  correlation,  here  again  there  is  difficulty  in  obtain¬ 
ing  the  accurate  relation  between  A  and  the  present  scale  of  turbulent 
vorticity  A. 

Besides  dealing  with  the  heat  conductivity  problem  in  turbulent 
flows,  the  paper  also  gives  applications  for  heated  wakes.  The  theo¬ 
retical  result  for  the  two-dimensional  heated  wake  agrees  with  the 
existing  experimental  measurement.’ 


Experimental  Verification  of  Several  Creep  Theories 
by  Creep  Tests  on  No.  45  Steel 

..  Utilizing  creep  curves  obtained  for  No.  45  steel  under  500°  C  and 
different  stresses,  several  typical  creep  theories,  i.e.-,  strain  harden¬ 
ing  theory — /(«■»</,«»)  =  0,  time  hardening  theory — /(*,<*,/)  =  «,  and 
plastic  hereditary  theory--  <iR(*)  =  <*(/)  +  —  r)*(f)</r,  are  verified.  • 

Results  are  analyzed  and  simple  methods  for  checking  these  theories-  are 
given. 
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Free  Edge  Rectangular  Plates  on  Elastic  Foundations 

This  paper  is  a  generalization  of  an  article  by  H.  Hapel  entitled 
"fiber  das  Gleichgevicht  von  elastischen -Flatten  unter  einer  Einzellast" 
(Mathematische  Zeitschrift,  6,  1920) .  It  differs  from  the  article  in 
four  respects: 

(1)  Both  symmetrical  and  asymmetrical  loadings  are  discussed  in 
detail. 

(2)  The  origin  of  the  coordinates  is  placed  at  the  corner  of  the 
plate  instead  of  at  the  center.  In  this  way  the  functions  U  and  V  and 
their,  integrals  are  considerably  simplified. 

(3)  Both  isotropic  and  orthotropic  plates  are  treated. 

(4)  Both  the  Ritz  and  the  Galerkin  method  of  variation  are  used 
to  solve  the  problem. 

■  In  the  paper  the  variation  equations  are  written  in  standard  form 
and  the  values  of  the  functions  U  and  V  and  their  second  derivatives  at 
ll  X  11  points  are  tabulated.  Once  the  dimensions  a  and  b,  the  rigidity 
D  ttf  the  plate,  the  modulus  k  of  the  foundation,  and  the  concentrated 
ib&ds  R  at*e  known,  the  reader  may  calculate  the  deflections  of  and  the 
moliieiits  ih  the  plate  by  direct  application  of  these  equations. 


Effect  of  Earthquakes  on  Dams  (l) 


In  this  paper,  the  results  of  a  model  test  of  earthquake  effect 
on  dams  are  presented. 

The  experimentation  consists  of  measuring  the  strains  and  displace¬ 
ments  at  various  points  on  the  surface  of  the  model  when  it  is  suddenly 
set  into  vibration  by  a  horizontal  impulse  applied  to  its  base.  The 
model  is  fixed  to  a  rigid  table  suspended  on  cables,  and  the  impulse  is 
calibrated  so  that  a  known  amount'  of  momentum  is  transferred  to  the 
table  by  impact  of  the  pendulum.  The  duration  of  impact  is  made  much 
shorter  than  the  second  natural  period  of  vibration  of  the  model.  The 
swing  of  the  table  after  impact,  being  of  very  low  frequency,  does  not 
cause,  any  appreciable  additional  distortion  to  the  model  and  can  be 
neglected;  . 

Insofar. as  the  model  may  be  regarded  as  a  linear  elastic  system, 
if  we  denote  the  stress,  (or •  displacement)  measured  at  any  point  by  K(t), 
then  for. arbitrary  ground  acceleration  2(0,  such  as  produced  by  an 
earthquake,  the  stress  (or  displacement)  at  the  same  point  may  be  eval¬ 
uated  by  means  of  Duhamel's  integral,  i.e., 

«K<)  =  K(j  —  r)  2(r) dx. 

The  model  is  made  of  rubber.  Bonded  resistance  wire  gauges  are 
used  for  the  measurement  of  strain.  The  strain  gauges,  because  of  their 
high  rigidity  in  comparison  with  that  of  rubber,  alter  in  a  significant 
way  the  local  strain  distribution  of  the  model,  so  .that  the  apparent  „ 
strain  readings  taken  by  the  ordinary  method  do  not  give  the  actual 
strains  which  would  be  produced  when  no  strain  gauges  were  used. 


However,  it  is  found  by  our  experiment  that  a  definite  ratio  exists 
between  the  apparent  and  actual  strains. 

The  results  of  the  experiment  are  given  in  Figures  3  to  8  (original 
text) . 

K(t)  can  be  resolved  into  various  components  of  free  vibration,  each 
component  being  of  the  type  sin  v  1  —  n)  mt  •  Conversion  of  the  result 

of  the  model  experiment  to  the  prototype  is  direct  provided  one  makes  cer 
tain  modifications  to  account  for  the  difference  in  damping'  coefficients 
and  neglects  the  difference  in  Poisson's  ratios.  Making  use  of  this  fact 
the  stresses  are  computed  by  means  of  Duhamel's  integrals,  for  which  an 
actual  acceleration  seismogram  is  used.  Similarly,  use  is  made  of  pub¬ 
lished  strong  earthquake  spectra  to  find  the  upper  bounds  of  stresses 
which  may  possibly  occur  during  strong  earthquakes. 

Efficiency  of  Water  Wheels  in  Hydrodynamics 

Among  the  various  small  types  of  low  water-head  power  systems  or 
hydroelectric  power  stations,  the  one  operated  by  water  wheels  is  the 
most  suitable  for  the  present  rural  areas  of  China.  The  use  of  water 
wheels  in  these  facilities  is  more  efficient  than  any  type  of  water  tur¬ 
bine,  especially  when  the  available  water-head  is  less  than  two  meters. 

In  this  paper  the  head  losses  in  the  "middle-push"  type  water 
wheel  are  analyzed,  their  magnitudes  estimated  and  the  similarity 'law 
derived.  It  is  pointed  out  that  the  kinetic  energy  flowing  out  from  the 
water  wheel  can  be  recovered  by  a  well-designed  open  channel. 
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ABSTRACT 


Design  for  the  Runner  Vanes  of  Francis  Water  Turbines  (l) 


Under  Vallander's  assumption  the  stream  surfaces  in  a  Francis  water 
turbine  remain  the  same  with  or  without-  the  runner  vanes.  The  equation 
governing  the  flow  in  a  thin  sheet  between  two  adjacent  stream  surfaces 

is 

■  i 

du1  dz2  A  dz  dz 


or 

d*tft  ,  d*</> _ 1  dh~  dip  __  0 

du  dz2  A  dz  dz 


where  <#*  and  4>  denote  the  velocity  potential  and  the  stream  function, 

••  and  z  are  the  transformed  coordinates  and  A  is  the  thickness  of  the 
thin  stream  surface  sheet  under  consideration. 

By  the  transformations  <f>  =  \/  A  <P  and  </>  =  — j=- &  >  the  pair  of  equa- 

v  A 

tions  is  reduced  to 


V2<£  =s  o  and  V2'f  =  o , 


are  small  by  comparison  with  the  square  of 


the  number  of  runner  vanes. 

The  flow  around  the  vanes  in  a  thin  stream  surface  sheet  was  also 
treated  by  Vallander  and  Kiselev;  however,  the  relation  between  the  flows 
in  two  adjacent  stream  surface  sheets  remains  as  yet  unknown.  In  this 
paper  the  same  relation  is  adopted  to  connect  the  flows' in  adjacent 
sheets  as  in  Bauersfeld-Voznesensky1 s  two  dimensional  theory,  that  is," 
the  circulation  along  the  contour  between  any  two  meridianal  planes  and 


IS 


outside  of  the  runner  vane  is  constant. 

The  solution  is  based  on  the  method  of  singularities.  At  the  end 
of  the  paper  the  vortex  distribution  function  is  discussed.  It  is  found 
that  the  addition  of  a  term  corresponding  to  the  flow  about  a  flat  plate 
results  in  a  definite  angle  of  attack  and  thereby  increases  load  near 
the  leading  edge  of  a  vane,  which  may  be  desirable  in  certain  casos.  But 
the  distribution  function  needs  to  be  further  improved  for  a  good  design. 

It  is  worthy  of  note  that  the  present  solution  is  also  applicable 
to  the  case  where  A  is  not  constant  and  small  a9  assumed  in  Kiselev's 
mothod,  and  that  it  would  avoid  the  longthy  computation  of  Vollnndor,i 
method,  although  Vallander's  solution  is  more  accurate. 


